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The E895 Collaboration at the Brookhaven ACS has performed a systematic inves- 
tigation of Au+Au collisions at 2-8 ACeV, using a large-acceptance Time Projection 
Chamber. In addition to extensive measurements of particle flow, spectra, two-particle 
interferometry, and strangeness production, we have performed novel hybrid analyses, in- 
cluding azimuthally-sensitive pion HBT, extraction of the six-dimensional pion phasespace 
density, and a first measurement of the A-p correlation function. 



1. Introduction 

Modern relativistic heavy ion physics focuses largely on the creation and study of the 
Quark Cluon Plasma (QCP), a state of matter in which partonic degrees of freedom are 
relevant over large length scales; potential signatures of QCP formation (see, e.g. [0]) 
are studied in collisions at ACS, SPS, and now RHIC. However, it is now well-accepted 
that only through a simultaneous study of several observables will a clear understanding 
of the collision history emerge. Furthermore, the most convincing evidence of "new" 
physics stems from a break in well-understood (in terms of hadronic models) systematics 
in observables which have been measured over a large range of conditions. 
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Using the EOS Time Projection Chamber (TPC) at the Brookhaven AGS, we have 
studied Au+Au colhsions at bombarding energies of 2, 4, 6, and 8 AGeV. Models sug- 
gest that colhsions at these energies reach the highest baryon density 0-0, making such 
systems ideal to stringently test our understanding of nuclear matter under extreme con- 
ditions, as well as to search for possible QGP formation in the high-density limit l^,^- 

We have measured excitation functions in four primary areas: (1) proton elliptic and 
directed flow (2) intensity interferometry [^-[1^, (3) charged particle spectra [p!3| , p!4 



and (4) strange particle production |jT5[. In addition to improving our understanding 
of hot hadronic systems under extreme conditions, these systematics provide a valuable 
baseline for measurements at higher energy. 

Going further, we have combined the analyses in these areas to achieve novel insights 
into the collision dynamics, and serve as an example of how to optimize physics informa- 



tion extraction for data at RHIC. In [^^-|T^, strange particle reconstruction is correlated 
with proton flow to extract information on in-medium potentials. In this report, we focus 
primarily on combined analyses involving two-particle correlations. 

2. Directed and Elliptic Flow 

Anisotropic flow refers to the momentum distribution of particles in azimuth about the 
beam axis. The reference azimuthal angle is the direction of the impact parameter vector 
b, which must be determined on an event-by-event basis; see for details. 

Two types of anisotropic flow are typically discussed [0. (Positive) sideward flow 



refers to the preferential in-plane emission of particles with y > ycm in the direction of b, 
and those with y < y^m opposite b, while elliptical flow analyses focus on the difference 
between in- and out-of-plane emission. In a simple picture in which the emission pattern 
can be described by an ellipsoid in momentum space (in general, the situation is more 
complicated [0), then, sideward flow analyses measure the tilt angle of the momentum- 



space ellipsoid, while elliptical flow probes the ratio between its semi-minor axes. 
E895 analyses of proton directed and elliptic flow have already been reported 



There, we found that AGS energies represent a transition region, in which the elliptic flow 
changes sign from out-of-plane (the "squeeze-out" signal observed at Bevalac/SIS energies) 
to in-plane elliptic flow seen at SPS energies, and directed flow reaches a maximum and 
falls off quickly with beam energy. Pion directed flow is oriented in the direction opposite 
proton flow- i.e. the vr momentum-space ellipse tilts away from b for > while 
elliptic flow changes sign at Ebeam ~4 AGeV, similar to the protons pT| . 

Dynamic transport models reproduce the gross trends in flow measured by E895, but 
the details are strongly sensitive to the underlying physics (e.g. the Equation of State) 
assumed hence, these excitation functions provide a very stringent test of our un- 

derstanding of the dynamics of heavy ion collisions. In particular, no single model is 
currently able to reproduce the detailed flow excitation function at AGS energies. 

Clearly, fresh information on the physics of anisotropic flow is needed. It is commonly 
assumed that these flow signals result from space-momentum correlations generated in 
the collision. By measuring azimuthally-sensitive momentum distributions, therefore, we 
look at only half the story. In the next section, we describe how two-pion correlations 
provide a more complete picture. 
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3. Pion Intensity Interferometry 

Intensity interferometry (a.k.a. HBT) techniques have been used extensively to probe 
the space-time structure of heavy ion colhsions 



22|,p3|). Until recently, most 



l^see, e.g. 

HBT analyses implicitly assumed a cylindrically symmetric emission region, strictly valid 
only for |b| = collisions. In this standard framework, E895 vr" correlation systematics as 
a function of beam energy, rriT, centrality, and rapidity, have already been reported P,p!2|. 



In addition to the (momentum-space) anisotropic flow studies discussed above, theo- 
retical studies p4|-[29[] suggest that coordinate- space anisotropics of the emitting source 
are equally interesting. Recently, a formalism was proposed to use azimuthally-sensitive 
intensity interferometry, correlated event-by-event with the reaction plane, to extract this 
information experimentally We have performed the first experimental measure- 

ment of the full coordinate-space anisotropics of a hot nuclear source, for pions emitted 
from Au-|-Au collisions at 2-6 AGeV. Full details may be found in [|T0[] . 

71" correlation functions were constructed in the Bertsch-Pratt ("out-side-long") de- 
composition as measured in the Au+Au cm. frame |30[| , 
with respect to the reaction plane = Z(K_|_,b), where K 
momentum of the pair perpendicular to the beam. 

The correlation functions for each of eight 45-degree-wide bins, are fit with the stan- 
dard Gaussian parameterization 

(1) 



with cuts on the pair angle 
± = (Pi + P2)± is the total 



C(q,0) = i + A(0)exp[- m3Rl{4>) 



i,j=o,s,l 

In contrast to 0-integrated analyses, all six HBT radii are relevant here P^P7| , PS| . 

For 4 AGeV collisions, the seven fit parameters R^j and A are plotted as a function of 
(j) in Figure |l|. The 0-independence of A indicates that the fraction of vr" from long-lived 
resonances is independent of emission angle with respect to the reaction plane, consistent 
with the relatively weak anisotropic fiow at these energies 0, as well as the randomizing 
effect on the pion from decay of the parent. Rq and Rg show significant equal and opposite 
2"'^-order oscillations in 0, while R^i and R^i display equal-magnitude l'^*-order oscillations. 

From the azimuthal dependence of the HBT radii, one may extract the geometry of 
the source in the impact-parameter-fixed coordinate system 
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In particular, at 

midrapidity and low px, R'ij{4') are sensitive to the homogeneity lengths in x, y, z, and 
time, as well as the tilt of the source in the reaction plane. Thus, we determine the 
shape and orientation of the freeze-out ellipsoid in coordinate space. The lines in Figure |I| 
represent a simultaneous fit to all six Rff-{(j)), using the formalism of p8 
following homogeneity lengths: ax 



resulting in the 
a, = 5.2 ± 0.1 



= 5.2 ± 0.1 fm, 
: 37° ±4°. Finite reaction-plane 
are accounted for 
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: 4.0 ± 0.1 fm, ay 

fm, at = 3.7 ± 0.5 fm/c; the extracted tilt angle is Og = 
resolution effects, which wash out oscillations in i?^^(0). 

Similar shapes and orientations were found for collisions at 2 and 6 AGeV. At all 
energies, our results indicate a pion freezeout distribution as an ellipsoid whose major 
axis in the reaction plane is tilted with respect to the beam in the positive direction (i.e. 
in the direction of b), and whose transverse axis perpendicular to the reaction plane is 
longer than the axis in the reaction plane. The extension in the temporal direction is 
consistent with observations at high energy 



,31 



The coordinate-space analog of elliptical flow, ay/a^, is shown as a function of beam 
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Figure 1. Azimuthal dependence of Bertsch- 
Pratt HBT parameters for midrapidity 
Tx~ from peripheral Au+Au collisions at 
4 AGeV. Lines represent a simultaneous fit 
to extract the underlying geometry. 



Figure 2. Excitation functions of 
the spatial analogs of elliptic (top) and 
directed flow (bottom) for data (filled 
stars) and RQMD (open symbols) with 
and without meanfield. 



energy in Figure 0. In contrast to the momentum space trends, there is no transition 
from out-of-plane to in-plane shape. Indeed, the observed transverse shape is reminiscent 
of the entrance channel geometry. In the simplest picture, pions are emitted from the 
overlap region of the two spherical Au nuclei. For impact parameter 6 = 4 — 8 fm, the 
ratio of out-of-plane to in-plane RMS for this overlap region varies from 1.3-1.9 (indicated 
by dashed lines in Figure consistent with the measured value ~ 1.3. 

The tilt angle, shown in the lower panel of Figure H, is the coordinate-space analog 
of directed flow, and is a geometric feature of the collision dynamics Studies with 
the RQMD model indicate that most of the low-pr pions arise from A decay, so it is not 
surprising that the tilted freeze-out distribution resembles the baryonic distribution in 
coordinate space calculated in transport codes . 

Experimental access to this level of geometric detail on the freezeout distribution is 
unprecedented, and represents an exciting new opportunity to study the dynamical re- 
sponse of hot nuclear matter to compression. An azimuthally-sensitive interferometry 
analysis [|10| on pions generated by the RQMD model displays considerable sensi- 
tivity on the dynamical effect of the nuclear meanfield. Although at these energies the 
model underpredicts the spatial and overpredicts the temporal scale (noted already for 
central collisions |^), qualitatively, the transverse shape and large positive tilt angles are 
reproduced remarkably well; transverse asymmetries and spatial tilt angles are shown in 
Figure]^. Since RQMD better describes proton directed flow (momentum-space tilt) when 
the meanfield is included in the calculation it is interesting to note that the spatial 
tilt angles are better described when the meanfield is off. Clearly, probing asymmetries in 
both momentum and coordinate space results in a more stringent test of the dynamics. 
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These coordinate-space anisotropics represent new and independent information. Mom- 
entum-space tilt angles (flow angles |^) at these energies are only a few degrees [Q, 



and, as mentioned above, for pions point in the negative direction opposite the 
coordinate-space tilt 6s- Experimental information on the interplay between coordinate 
and momentum space anisotropics should help resolve theoretical issues, such as the 
coexistence of flow and antiflow components |26[ | . 

4. Six-dimensional Pion Phasespace Density 

We have measured an excitation function of charged particle spectra over a wide region 
of phasespace. Figure |^ shows rriT distributions for tt" from central Au+Au collisions [r^|. 



Assuming a thermal picture, these spectra may be combined with those of heavier parti- 
cles, to extract freezeout temperatures and radial flow values 



As a test of this thermal assumption, it has been observed that for a thermally equi- 



librated pion source, both the shape and the amplitude of the six- dimensional phasespace 
occupation density as a function of spatial and momentum coordinates, f{x,p), are deter- 
mined by the temperature. Experimentally, we may probe the spatially-averaged phases- 
pace density by combining spectral and interferometry information, according to p^p5| 



1 d^N 



V Ro{pT,y) ■ Rs{PT,y) ■ Ri{PT,y) 

(For noncentral collisions or away from midrapidity, "cross-term" HBT radii also enter 
the denominator.) The factor from the HBT analysis is included in Eq. ^ to account 
for pions from long-lived decays [^,^. 

For heavy ion collisions at maximum AGS and SPS energies, a remarkable "universal" 
phasespace density {f){pT,y) has been observed ^Sf, independent of collision energy or 



system, though multiplicities vary by more than an order of magnitude. This suggests 
that vr — vr cross-sections dominate the freeze-out density of high-energy collisions, and it 
is natural to look for a break-down of the universal behavior at lower energy. 

Figure ^ shows the n~ average phasespace density {f){pT, 1/ = 0) determined according 
to Equation ^ for central Au+Au collisions at 2-8 AGeV. Clearly, the universal behavior 
has broken down; (/) grows monotonically with collision energy at the AGS. As at higher 



energies, (/) <^ 1, indicating that multiparticle symmetrization effects are weak p3| , |35 
Lines indicate expectations for "pure thermal" (no flow) sources with T ~ 100 MeV; The 
extracted values appear roughly consistent with the shape and magnitude of a thermal 
source, with a temperature increasing with Ebeam- In detail, however, the observed {f){pT) 
is flatter than the "pure thermal" expectation, consistent with a high-p^ enhancement 
due to radial flow [p5|; preliminary fits (not shown) with a parameterization that includes 



radial flow |36] reproduce the distributions well, and indicate that both temperature and 
flow increase with beam energy. 

Also shown on Figure ^ are values from Pb+Pb collisions at 158 AGeV ||3^ and pre- 



liminary data from STAR (Au(130 AGeV)+Au(130 AGeV)) 0, demonstrating the per- 
sistence of the universal behavior to RHIC energies. While the phasespace densities 
achieved at SPS and RHIC are higher than those at AGS at higher pt, presumably due 
to the greater explosive flow generated, it is remarkable that the "universal" value at low 
Pt is reached already for collisions below maximum AGS energy. 
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Figure 3. Preliminary acceptance-corrected vr" spectra from central Au+Au collisions at 
2-8 AGeV, in 0.1-unit-wide rapidity slices. The m-r spectrum at midrapidity is shown 
by the open boxes; spectra for other rapidity slices are scaled by factors of ten. Lines 
represent a two-component thermal fit. 

5. A-p Correlations 

In the baryon-rich system generated at the AGS, it is interesting to study the baryon 
freezeout geometry. Two-proton correlations may be used to extract this information 



E895 proton interferometry results indicate a proton freezeout volume independent of 
beam energy |18]. It has been suggested that A — p correlation measurements could offer 
more sensitive determinations of baryon source sizes than pp interferometry, particularly 
for large sources. These measurements may also shed light on the A-p interaction, which 
is not fully understood PB| . 

For semi-central {b < 7 fm) Au-|-Au collisions, ~ 120, 000 As were identified in the 
TPC through their decay topology into p + tt~ , using a neural network |16,39]. The A- 
p correlation function, as a function of the relative momentum in the pair rest frame, 
is shown in Figure ^. A significant enhancement is seen at low relative momentum, 
consistent with an attractive A-p final state interaction. 

Also shown are expectations for the correlation function for various-sized static sources, 
using a parameterized final state potential . The shape of the experimental correlation 
function does not match these expectations for any source size. However, since the effects 
of the reconstructed A purity (~ 80%) and of A feed-down from the electromagnetic decay 
of S°s (estimated at 25 — 30%) have not yet been accounted for, it is premature to draw 
conclusions regarding source size or A-p potential. 
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Figure 5. Preliminary A-p correlation 
function for Au+Au collisions at 6 AGeV 
(a) and for combined data from collisions 
at 2, 4, and 6 AGeV (b). Also shown are 
expectations from a static source 
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Figure 4. Preliminary phasespace 

density for midrapidity tt~ from central 
Au(Pb)+Au(Pb) collisions, measured by 
E895, NA49 [H, and STAR 13- Uncertain- 
ties, not shown, are 10-20%. 

6. Summary 

The E895 Collaboration has performed extensive studies of spectra, interferometry, 
strangeness, and flow at AGS energies, providing valuable energy systematics useful for 
identification of "new" physics at SPS and RHIC. Further, we have combined these anal- 
yses to extract new physics information. The azimuthal dependence of pion HBT radii 
reveals the anisotropic transverse shape of the emitting source, and its tilted orientation 
in the reaction plane, imaging proton and pion flow in coordinate space for the first time. 
The 6-dimensional vr" phasespace occupation density is consistent with emission from a 
flowing, thermal source; {f)^- grows with beam energy, approaching the "universal" limit 
observed at SPS and RHIC energies. First measurements of the A-p correlation function 
display an enhancement at low relative momentum; more detailed study of this correla- 
tion should provide valuable information on the baryon freezeout geometry and the A-p 
potential. 
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